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ABSTRACT 

We use -88 arcmin^ of deep (>26.5 mag at 5cr) NICMOS data over the two GOODS fields and 
the HDF South to conduct a search for bright z>7 galaxy candidates. This search takes advantage 
of an efficient preselection over 58 arcmin^ of NICMOS i^ieo-band data where only plausible z ^ 1 
candidates are followed up with NICMOS Jno-band observations. —248 arcmin^ of deep ground- 
based near-infrared data (> 25.5 mag, 5<j) is also considered in the search. In total, we report 15 
Z85o-dropout candidates over this area - 7 of which are new to these search fields. Two possible 2: — 9 
Jiio-dropout candidates are also found, but seem unlikely to correspond to z — 9 galaxies (given 
the estimated contamination levels). The present z — 9 search is used to set upper limits on the 
prevalence of such sources. Rigorous testing is undertaken to establish the level of contamination 
of our selections by photometric scatter, low mass stars, supernovae (SNe), and spurious sources. 
The estimated contamination rate of our z ^ 1 selection is —24%. Through careful simulations, the 
effective volume available to our 2: > 7 selections is estimated and used to establish constraints on 
the volume density of luminous {L*^^,^ ——21 mag) galaxies from these searches. We find that 
the volume density of luminous star-forming galaxies at z — 7 is 13^5 x lower than at z — 4 and 
>25x lower (la) at z — 9 than at z — 4. This is the most stringent constraint yet available on the 
volume density of > ^^=3 galaxies at z — 9. The present wide- area, multi- field search limits cosmic 
variance to <20%. The evolution we find at the bright end of the UV LF is similar to that found from 
recent Subaru Suprime-Cam, HAWK-I or ERS WFC3/IR searches. The present paper also includes a 
complete summary of our final 2: — 7 zgso-dropout sample (18 candidates) identified from all NICMOS 
observations to date (over the two GOODS fields, the HUDF, galaxy clusters). 

Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

The recent WFC3/IR camera on the Hubble Space 
Telescope (Kimble et al. 2006) has completely revolution- 
ized our ability to search for galaxies at z > 7 due to its 
extraordinary imaging capabilities in the near-infrared - 
allowing for large areas to be surveyed to great depths. 
Already some 40 credible z— 7-8 galaxy candidates have 
been identified in the ffist hundred orbits of observations 
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(Oesch et al. 2010a; Bouwens et al. 2010a, b; McLure et 
al. 2010; Bunker et al. 2010; Yan et al. 2010; Finkel- 
stein et al. 2010; Wilkins et al. 2010a,b). This com- 
pares with —15 credible candidates reported thus far 
from deep, wide-area ground-based observations (Ouchi 
et al. 2009; Castehano et al. 2010; Hickey et al. 2010) and 
-12 thusfar with NICMOS (e.g., Bouwens et al. 2008, 
2009a; Bradley et al. 2008; Oesch et al. 2009; Zheng et 
al. 2009; - 2 from Richard et al. 2008). Whereas -100 
orbits of NICMOS observations were required to obtain 
1 2: — 7 credible galaxy candidate (see e.g. §4 of Bouwens 
et al. 2009a), only -2.5 orbits of WFC3/IR observations 
are required to find a similar 2^-7 candidate. 

Despite these significant advances in our observational 
capabilities with WFC3/IR to reach deep and identify 
large numbers of faint z^l galaxies, a full characteriza- 
tion of the galaxy population at z — 7 requires that we 
identify large numbers of galaxies at both high and low 
luminosities. All but —6 galaxies in early selections of 
z— 7-8 galaxies from early WFC3/IR observations over 
the ERS/HUDF09 fields have magnitudes faintward of 
26.5 mag (e.g., Oesch et al. 2010a; Wilkins et al. 2010a; 
Bouwens et al. 2010c). As such, it is somewhat challeng- 
ing to characterize the properties of relatively luminous 
galaxies at z > 6.5, and some expansion of the number 
of sources known brightward of 26.5 mag would be ben- 
eficial. Such samples are particularly valuable over fields 
such as GOODS (Giavalisco et al. 2004) where other valu- 
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able multiwavelength data exist like deep IRAC (Dickin- 
son & GOODS Team 2004) or Chandra coverage (Brandt 
et al. 2001; Rosati et al. 2002). 

Fortunately, for the selection of luminous z^7-8 galax- 
ies, some '^88 arcmin^ of deep (>26.5 mag, 5cr), wide- 
area NICMOS observations exist over the two GOODS 
fields and the HDF South. While ~23 arcmin^ of those 
observations have already been used to identify z^7-S 
galaxies (Bouwens et al. 2008; Oesch et al. 2009), >60 
arcmin^ of those data have yet to be used. Most of these 
data are associated with the GOODS NICMOS survey 
(Conselice et al. 2010) or are NICMOS parallels asso- 
ciated with other HST programs. In total, these NIC- 
MOS observations cover 2x as much area as available in 
the WFC3/IR observations that made up the Early Re- 
lease Science Program (PI O'Connell: GO 11359). Mean- 
while, these observations cover comparable area to that 
available in the wide-area HAWK-I F-band observations 
over the CDF-South (Castellano et al. 2010; Hickey et 
al. 2010) and less area than the Subaru Suprime-Cam 
(Ouchi et al. 2009) observations over the HDF-North, 
but are deeper on average (by '^0.2 and ~0.7 mag, re- 
spectively). 

In the present work, we use these wide-area NICMOS 
observations to identify a small sample of luminous star- 
forming galaxies at z ~ 7, adding significantly to that 
known, and performing a search at z ~ 9. Some of the 
present z ^ 7 sample has already been used by Gonza- 
lez et al. (2010) and Labbe et al. (2010b) to do stellar 
population modelling of luminous z ^ 7 galaxies and to 
extend measures of the stellar mass density and specific 
star formation rate to z ^ 7. Here we describe the se- 
lection of those z > 7 candidates in detail, discuss the 
properties and layout of the NICMOS observations, and 
summarize the properties of the sample. We also esti- 
mate the contamination rates and the selection volume 
for this sample. Finally, we use this search to derive a 
constraint on the volume density of L^^^ (or ~— 21 mag) 
galaxies at z ~ 7. We will also incorporate results from 
the ~248 arcmin^ Bouwens et al. 2008 search for 2: ~ 7 
galaxies in deep ground-based data over GOODS. 

The structure of this paper is as follows. In §2, we 
summarize the observational data. In §3, we describe 
our search results for z > 7 galaxies. In §4, we use 
the observational search results to derive a constraint 
on the bright end of the z ^ 7 and z ^ 9 UV LFs. 
Finally, we provide a summary (§5). Throughout this 
work, we often denote luminosities in terms of the char- 
acteristic luminosity I/J=3 at 2; ^ 3 (Steidel et al. 1999), 
i.e., Muv,AB = —21.07. Where necessary, we assume 
Qo = 0.3', = 0.7, Ho = 70km/s/Mpc. Although 
these parameters are slightly different from those de- 
termined from the WMAP seven-year results (Komatsu 
et al. 2010), they allow for convenient comparison with 
other recent results expressed in a similar manner. We 
express all magnitudes in the AB system (Oke & Gunn 
1983). 

2. OBSERVATIONAL DATA 

The present search for galaxy candidates at z > 7 
makes use of --88 arcmin^ of deep NICMOS + ACS ob- 
servations over the two GOODS fields + HDF South, as 
weh as --248 arcmin^ of deep MOIRCS+ISAAC ground- 
based data. 



Table [T] provides a convenient summary of the observa- 
tional properties of each of these search fields, while Fig- 
ure [T] shows their position within the two GOODS fields. 
A significant part of these observations (e.g., those over 
the HDF North [Wilhams et al. 1996], HUDF [Beckwith 
et al. 2006], NICMOS parallels to HUDF) have already 
been used to search for 2; > 7 galaxies (e.g., Bouwens et 
al. 2008; Oesch et al. 2009), and so we will not repeat a 
detailed description of those observations here. 

Instead we focus on the NICMOS observations taken 
over the past three years not yet utilized for z > 7 
searches. These observations include 317 orbits of NIC- 
MOS data from five different HST programs (G09723, 
GO10403, GO11082, G011144, G011192) over >60 
arcmin^. 

The foundation of our search for new z ^ 7 candi- 
dates is provided by the NICMOS observations of the 
first three programs (G09723, GO10403, and GO11082). 
These programs involve 236 orbits of i^ieo-band NIC- 
MOS imaging observations, and cover --58 arcmin^ of 
the —65 arcmin^ of new search areaO 180 of these or- 
bits came from the GOODS NICMOS Survey (Conselice 
et al. 2010), with 60 separate 3-orbit NIC3 pointings. 40 
of the orbits were obtained in parallel with ACS SBC 
far-UV observations over the HUDF (Siana et al. 2007; 
GO10403: PI Teplitz), and 16 of the orbits came from 
a NICMOS program over the HDF-South WFPC2 field 
(Labbe et al. 2003; G09723: PI Franx), with 8 2-orbit 
NIC3 pointings. We have indicated the position of those 
NIC3 pointings that lie within the GOODS fields on Fig- 
ure [T] with light orange squares. These observations were 
reduced with our NICMOS pipeline "nicred.py" (Magee 
et al. 2007) and reach 26.7 and 26.9 mag (5cr, 0.6''- 
diameter apertures) depending upon whether the i^ieo- 
band integrations were 2 or 3 orbits, respectively, in du- 
ration. The FWHM for the iJieo-band PSF is -0.37''. 

The first two of these programs (GO 10403 and 
GO 11082) had deep optical coverage from the GOODS 
ACS program, with corresponding depths reaching 27.9, 
28.1, 27.4, and 27.3 in the ^435, V^eoe, ^775, and Zs50 
bands, respectively, in 0.6"-diameter apertures. Our re- 
ductions of the ACS GOODS observations are described 
in Bouwens et al. (2007) and take advantage of essen- 
tially all ACS data ever taken over the GOODS fields. 
They are similar to the GOODS v2.0 reductions, but 
reach —0.1-0.3 mag deeper in the 2:85o-band due to the 
inclusion of the SNe follow-up data (Riess et al. 2007). 
The third program here (G09723) had deep WFPC2 ob- 
servations from the HDF South program (Williams et al. 
2000), and 2 orbit ^775/2 orbit zs^o ACS WFC observa- 
tions as part of the ACS GTO program (PI Ford: GTO 
9301). 

The IRAC reductions we use for our selection of z ^7 
galaxy candidates are those from Dickinson & GOODS 
Team (2004). They reach to 27.4 mag, 26.8 mag, 25.4 
mag, and 25.3 mag in the 3.6/im, 4.5/im, 5.8/im, and 
8.0/im bands, respectively, at la (2") in most regions 
(where 23 hour integrations are considered) and 0.4 mag 
deeper (where the integration time is 46 hours). 

Deep Jiio-band observations became available over 
some fraction of these fields, as a result of our 60-orbit 

We note that a small fraction (~3-4 arcmin-^) of this search 
area had been previously considered (Henry et al. 2009). 
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Fig. 1. — Deep, wide-area near-IR data available to search for z^7 galaxies over the CDF-South GOODS {left) and HDF-North GOODS (right) 
fields (see also Table fT]) . Most of the ultra-deep NICMOS and deep NICMOS observations were already presented in Bouwens et al. (2008) and 
included here on this figure (red and dark orange regions corresponding to regions with 5cr depths of > 28 mag and > 26.5 mag, respectively). 
What is new in the current analysis is the moderately deep (26.9 mag at 5a), wide-area (~58 arcmin^: >60 NIC3 pointings) NICMOS i^ieo-band 
data (light orange squares: Conselice et al. 2010). These data can be used to search for 2 > 7 galaxy candidates by identifying those sources that 
are not detected in the ACS optical bands, but have very red Z850 — Hieo colors and possess reasonably blue H — 5.8/im and 3.6/xm — 5.8/im colors. 
Unfortunately, these criteria (while very demanding) are not sufficient to place strong enough constraints on the redshift of the candidates, and so 
we also obtained deep (27.0 mag at 5a) NICMOS Jno-band imaging (magenta squares) over the best 2; > 7 candidates f Table [2] see §3.2) with 
GO program 11144 (PI: Bouwens). Also shown in dark orange and annotated with "Y" or "H" (for the G011192 H. Yan et al. [2010, in prep] 
or Henry et al. 2009 fields, respectively) are several additional NICMOS fields we used to search for z ^ 7 galaxies not considered by Bouwens 
et al. (2008). NICMOS fields previously considered by Bouwens et al. (2008) in searches for 2 > 7 galaxies are indicated as follows: "D" denotes 
the HDF-North Dickinson field (Dickinson 1998), "T" denotes the HDF-North Thompson field^(Thompson et al. 1999), "U" denotes the HUDF 
Thompson field (Thompson et al. 2005), "S" denotes the HUDF05 NICMOS field over the HUDF (Oesch et al. 2009), "1" denotes the first set of 
NICMOS parallels to the HUDF (NICP12: Oesch et al. 2009), "2" denotes the second set of NICMOS parallels to the HUDF (NICP34: Oesch et al. 
2009), and "G" denotes the GOODS Parallel NICMOS fields. The blue and dark blue regions correspond to regions with deep and very deep optical 
ACS coverage, respectively {5a depths of ^ 28 and > 29 mag). The position of the ground-based ISAAC+MOIRCS search areas is not shown here 
to minimize confusion (but is presented in Figure 1 of Bouwens et al. 2008). Also not included on this figure are the NICMOS search fields (Zirm 
et al. 2007) over the HDF South (Williams et al. 2000). The position of the Early Release Science WFC3/IR observations with the CDF South 
(not used here for 2; ~ 7 LF constraints) is indicated by the light shaded red region. A link to our NICMOS reductions over the GOODS fields is 
provided at http://firstgalaxies.0rg//astr0n0mers-area 



G011144 NICMOS program. The Jno-band observa- 
tions were '^2-3 orbits in duration, reaching depths of 
26.9-27.0 mag (5cr in 0.6''-diameter aperture) and tar- 
getted specific z ^7 candidates identified in the original 
ACS optical + NICMOS i^ieo band observations (see 
§3.2 for a description of the preselection). Only 14 such 
NIC3 fields were obtained and are indicated in Figure [T] 
with the magenta squares (though two of the NIC3 fields 
shown targetted specific z^6 candidates to allow for a 
measurement of the t/y-continuum slope /3: see Bouwens 
et al. 2009b). Again the reduction of those observations 
was performed with "nicred.py" (Magee et al. 2007), an 
alignment made to the ACS GOODS data, and then the 
data drizzled onto the ACS GOODS frame rebinned on a 
0.09'' X 0.09'' frame. The FWHM for the Jno-band PSF 
is -0.33". 

Two of our NIC3 fields chosen for Jno-band follow-up 
were allocated more than 3 orbits of time, to better ascer- 
tain the nature of the z >! candidates in those fields. In 
the first case (S4: 03:32:42.48, -27:42:48.8), additional 
follow-up observations were obtained with NICMOS. In 
the other case (S2: 03:32:06.75, -27:47:07.0), additional 
follow-up observations were obtained with WFC3/IR in 
the Ji25 band (1 orbit). 



While most of the search power of our program is pro- 
vided by the first three iJieo-band programs and the 
follow-up Jno-band program, we also include observa- 
tions from two other HST programs in our search for 
z > 7 galaxies. One of these programs is the H. Yan et al. 
(2010, in prep) GO 11192 program designed to follow up 
on bright zgso-dropout candidates identified in wide-area 
ground-based near-IR observations over the two GOODS 
fields (see the dark orange squares in Figure [T] annotated 
with "Y"). In that program, 4 orbit Jno + iJieo band 
observations (2 orbits in the Jno band and 2 orbits in 
the HiQo band) were obtained on 6 different bright z ^ 7 
candidates. These fields reached depths of 26.9 and 26.7 
mag in the Jno and iJieo bands, respectively, in 0.6"- 
diameter apertures. 

The other program that we utilized was GO 10872 
(Henry et al. 2009; Siana et al. 2010). NICMOS par- 
allel observations from that program cover —3.4 arcmin^ 
of area and lie within the two GOODS fields (see the dark 
orange squares in Figure [T] annotated with "H"). These 
observations reach to —26.9 mag in the Jno band and 
-26.7 mag in the iJieo band. While Henry et al. (2009) 
have already used these fields to search for candidate 
z > 7 galaxies - reporting none - we shall nevertheless 
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TABLE 1 

ACS + NICMOS + Ground-Based Imaging data used for our z and J 

DROPOUT SEARCHES.^ 







5cr Depth^ 






Name 


Area zsso 


-^110 ^160 


Ks 


Ref^ 



New NICMOS Fields 



GOODS NICMOS Survey 


44.6 


27.5 


d 


26.9 


^25^ 


[1 




Teplitz Parallels^ 


9.3 


27.5 


d 


26.9 


-25^ 


[2 




HDF-South 


4.3 


26.4 


26.3^ 


26.7 


25.7 


[3: 




Yan Survey 


3.1g 


27.5 


26.9 


26.7 




[4 




Henry Fields 


3.4g 


27.5 


26.9 


26.7 




[5] 




NICMOS Fields Already Considered 


by Bouwens et al. 


(2008) 




HDF-North Dickinson 


4.0 


27.8 


27.0 


27.0 


25.6 


[6,7] 


HDF-North Thompson 


0.8 


27.8 


28.0 


28.1 


25.6 


[8,7] 


HUDF Thompson 


5.8 


29.0 


27.6 


27.4 


26.0 


[9,10] 


HUDF Stiavelli 


0.7 


29.0 


28.1 


27.9 


26.0 


10,11 


HUDF-NICPARl 


1.3 


28.6 


28.6 


28.4 




11,12 


HUDF-NICPAR2 


1.3 


28.6 


28.6 


28.4 




11,12^ 


GOODS Parallels 


9.3 


27.5 


27.0 


26.9 


- 25^ 


[12] 


Fields with Ground-Based Data Already Considered 


in Bouwens et al. (2008) 




ISAAC v2.0 (CDFS) 


136 


27.5 


- 25.4^ 


- 24.8^ 


- 25^ 


[13] 


MOIRCS GTO-2 (HDFN) 


28 


27.5 


25.6 




25.6 


[7] 


MOIRCS GTO-1,3,4 (HDFN) 


84 


27.5 


24.2 




24.4 


[7] 



^ The layout of these search fields is illustrated in Figure [T] 

^ 5cr depths for ACS and NICMOS data given for a 0.6^^-diameter aperture and for a 
~ 1.0^^-diameter aperture for the ground-based Kg-band data. No correction has been 
made for the nominal light outside these apertures (for example, for a point source, 
the correction is typically ~0.2 mag) to keep the present estimates as empirical as 
possible. This is in contrast to the convention that we use in some previous work 
(e.g., Bouwens et al. 2008) where such corrections have been made. 
^ References: [1] Conselice et al. 2010, [2] Siana et al. (2007), [3] Labbe et al. (2003), 
Zirm et al. (2007) [4] H. Yan et al. (2010, in prep) [5] Henry et al. (2009), [6] Dickinson 
1998, [7] Kajisawa et al. 2006, Ouchi et al. 2007, [8] Thompson et al. 1999, [9] 
Thompson et al. (2005), [10] Labbe et al. (2006), [11] Oesch et al. 2007, [12] Bouwens 
& lUingworth (2006), Riess et al. (2007), Siana et al. (2007), [13] Retzlaff et al. 2010, 
Mannucci et al. 2007, Stanway et al. (2008). 

^ NICMOS Jiio-band observations, with 5(7 depths of 27.0 mag (0.6^^-diameter aper- 
ature) were acquired in those fields with 2; > 7 candidates. 

^ 40-orbit NICMOS ifieo-band observations taken in parallel with ACS SBC far-UV 
observations of the HUDF (Siana et al. 2007; GO10403: PI Teplitz) 
^ The J-band observations here are from the deep FIRES observations over the HDF- 
South with ISAAC (Labbe et al. 2003). 

g Not including overlap with the GOODS NICMOS Survey (Conselice et al. 2010) 
^ The depth of the near-IR data over the CDF-South varies by ~0. 2-0.4 mag depend- 
ing upon the observational conditions in which the ISAAC data were taken (Mannucci 
et al. 2006; Stanway et al. 2008; Retzlaff et al. 2010). 

be the square root of the image (Szalay et al. 1999). 
Here the square root of the image is simply the NIC- 
IMOS i^i6o-band imaging data. Colors were measured in 
a smaller scalable aperture using Kron-style (1980) pho- 
tometry, with a Kron factor of 1.2 to improve the S/N 
(typical aperture diameters for color measurements were 
0.4^'). To correct the small aperture flux measurements 
to total magnitudes, we computed corrections by com- 
paring the light in a larger scalable (Kron factor of 2.5) 
aperture on the image with that in the smaller aper- 
ture (Kron factor of 1.2). Figure 5 of Coe et al. (2006) 
provides a graphical description of a similar multi-stage 
procedure for measuring colors and total magnitudes. 
We also make a small correction for light outside the 
larger scalable apertures and on the wings of the PSF 
(typically ~0.15 mag). This correction is made based 
upon the encircled energy expected to lie outside this 
aperture (for stars). 



3.2. z^l Candidate Preselection 



incorporate these observations into our 2; > 7 search to 
be as comprehensive as possible. 

3. SELECTION OF Z > 7 CANDIDATES 

Here we describe the selection of z > 7 galaxies over 
those NICMOS fields not previously considered by our 
team for such searches (~65 arcmin^ in search area). We 
will not revisit the searches for z > 7 galaxies conducted 
by Bouwens et al. (2008) over the ~23 arcmin^ of area 
within the HUDF, HDF North, and CDF South or the 
~248 arcmin^ of deep ground-based data. Instead we 
will simply incorporate the Bouwens et al. (2008) search 
results in with our new search results. 

3.1. Catalogues 

We use the same procedure for doing object detection 
and photometry as we have used in previous work (e.g., 
Bouwens et al. 2007; Bouwens et al. 2008). Briefiy, we 
begin by PSF-matching all of our data to the PSF of 
our detection image (which here is the NICMOS i^ieo 
band). SExtractor (Bertin & Arnouts 1996) is then run 
in double- image mode, with detection image taken to 



5 



TABLE 2 

Z>7 CANDIDATES IN OUR WIDE-AREA NICMOS ifieo-BAND DATA PRESELECTED FOR Jhq-BAND FOLLOW-UP OBSERVATIONS.^ 



ID 


R.A. 


Dec 


NICMOS 
Follow- Up 
Pointing^ 


B-H 


V-H 


i-H 


z- H 


z- J 


J - H 


H 


(7{HY 








Z 


' 7 2;85o-dropout candidates 












GNS-zDl^ 


03:32:43.29 


-27:42:47.9 


S4 


>3.2 


>3.7 


2.7±0.7 


1.7±0.3 


1.3±0.4 


0.4±0.2 


25.8±0.2 


9 


GNS-ZD2 


03:32:32.03 


-27:45:37.2 


S5 


>2.6 


2.9±0.6 


>2.3 


>2.2 


1.6±1.0 


0.6±0.3 


26.2±0.3 


5 


GNS-zD3^'f 


03:32:06.10 


-27:46:37.3 


S2 


>3.1 


>3.4 


>2.9 


>2.5 


1.5±1.0 


0.8±0.5 


26.1±0.3 


5 


GNS-ZD4 


12:36:10.93 


62:09:15.6 


N3 


>2.6 


>3.0 


>2.4 


>2.3 


>1.6 


0.7±0.3 


25.8±0.3 


6 


GNS-zD5f 


12:36:44.68 


62:16:15.4 


N4 


>3.3 


>3.6 


>3.2 


>2.5 


1.7±0.8 


1.0±0.3 


25.4±0.2 


10 


GNS-zD6^'S 


03:32:22.66 


-27:43:00.6 


SI 


>3.2 


>3.4 


>2.6 


1.9±0.4 


1.6±0.5 


0.2±0.2 


25.4±0.2 


11 


GNS-zD7^'^ 


03:32:42.84 


-27:42:47.7 


S4 


>2.5 


>3.0 


>2.3 


1.9±0.5 


1.5±1.0 


0.4±0.3 


26.1±0.2 


5 








z 


~ 9 J-dropout candidates 












GNS- TD1 ^ 


n'^-'^9-1 77 

kjO.OAi. J.O. ( ( 




S8 




2.9±0.9 


2.5±0.9 


>1.9 


>0.5 


1.4d=0.5 


oc: 7-l-n o 

ZiO . I ZCU. Zi 


6 




1z:od:z5.47 


oz:14:ol.o 


N2 
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Probable Spurious J-dropout Candidates^ 












03:33:04.66 


— 27:52:27.0 


S3 




>3.4 


>3.2 


2.3±0.5 




>2 


26.5±0.3 


5 


GNS-Sp2 


03:33:04.14 


-27:52:57.4 


S3 




>3.4 


>3.0 


>2.3 




>2 


26.0±0.2 


5 


GNS-Sp3 


03:32:08.06 


-27:46:58.1 


S2 


>2.9 


>3.1 


>2.9 


>2.4 




>2.1 


27.0±0.3 


5 








Other Candidates Targetted 


in Preselection^ 










GNS-Ol 


03:33:02.18 


-27:53:40.4 


S6 


>2.0 


>2.2 


0.8±0.6 


1.7±0.8 


0.9±0.8 


0.7±0.5 


26.7±0.3 


5 


GNS-02 


03:32:27.81 


-27:54:48.0 


S7 


2.7±0.6 


>3.4 


>2.9 


1.4±0.4 


0.3±0.5 


1.1±0.4 


25.1±0.2 


7 


GNS-03 


12:36:15.21 


62:10:39.7 


N6 


>2.3 


>2.5 


>2.0 


1.4±0.5 


0.6±0.7 


0.8±0.6 


26.7±0.3 


5 


GNS-04 


12:36:11.26 


62:09:00.0 


N3 


2.1±0.6 


>2.8 


>2.2 


1.3±0.5 


0.3±0.6 


0.9±0.5 


26.3±0.3 


5 


GNS-05 


12:36:19.14 


62:15:23.4 


Nl 


2.8±0.9 


>3.0 


>2.6 


1.8±0.6 


0.6±0.7 


1.2±0.4 


25.5±0.2 


5 


GNS-06 


12:36:31.68 


62:06:45.8 


N5 


>2.7 


2.6±0.6 


>2.4 


>2.3 




>2.0 


26.2±0.3 


5 


GNS-07 


12:37:06.51 


62:11:49.0 


N7 


>2.8 


>3.0 


>2.6 


2.1±0.6 


<-0.2 


>1.9 


26.1±0.3 


5 


GNS-08 


12:37:03.01 


62:11:35.8 


N7 


>3.0 


>3.1 


>2.7 


1.2±0.3 


0.3±0.4 


0.9±0.3 


26.0±0.3 


6 



Lower limits are Icr. 

^ NIC3 pointing (from G011144) used for JnQ-hand follow-up observations on specific z > 7 candidates. Fourteen different NIC3 pointings were utilized (SI, 82, S3, 84, 
85, 86, 87, 88, Nl, N3, N4, N5, N6, N7). 

Significance (in cr) of the //j^gQ-band detection. The significance is measured in our smaller scalable apertures (Kron factor of 1.2). These apertures are smaller than 
those used for estimates of the total magnitudes (§3.1), and hence the associated significance levels are accordingly higher. 

Flux information is also available on these candidates from the WFC3/IR ER8 program (see Table [3j. 

^ WFC3/IR J]^25"t>^^d data is also available for GN8-zD3 as a result of our GO 11144 program. The measured J]^25"t>^^d magnitude for this candidate is 26.5ih0.2 mag. 

The Ji25 — -fi^ieO color (0.4±0.3) is somewhat red for a z ~ 7 galaxy. GN8-zD3 is thus more likely at low redshift than most candidates in our selection, 
f 

These z ~ 7 candidates are sufficiently red in their Jhq ~ -^^160 colors that there is an increased probability they could be low redshift contaminants. We have no 
evidence for this possibility however. 

3 Also identified in Castellano et al. (2010) and Hickey et al. (2010) 

^ This source has a measured ^850 ~ ^105 color ~ 0.0 it 1.0 and J±25 ~ ^160 ^ ^ ^-^ color in the high 8/N WFC3/IR observations (see Table[3j. This suggests this 
candidate may be a red z~l-2 galaxy and not a z ~ 7 star- forming galaxy. 

* GN8-JD1 is detected at Icr in the VgQg and ■1775 bands, not sufficient to rule it out as a z ~ 9 Jhq dropout candidate, but suggesting it may lie at z~l-2. 
3 We have no evidence that GN8-JD2 is detected at wavelengths other than l.G^m (irf-band). It is sufficiently close to another source that it is unclear if it is detected 
redward of 2^m from the IRAC data. This may suggest it is a transient source (SNe: see §3.4) or spurious (since it is close to the edge of the NIC3 field where it was 
found). However, we cannot rule out the possibility it corresponds to a z ~ 9 galaxy (but we consider it unlikely). 

^ These candidates were targetted for follow-up observations based upon ~5cr detections in the NICMOS /f^^gQ-band data. However, they are not detected in the NICMOS 
JWQ band observations. Lacking such detections, their reality is less secure. We suspect that they may be spurious (see §3.4). 

^ These candidates were targetted for J]^]^Q-band follow-up observations with NICMOS (because of their red zg5Q — -f^igo > 1-2 colors, apparent absence of flux in the 
optical, and moderately blue H — 5.8^m colors: see §3.1). However, they do not satisfy our z ~ 7 zg5Q-dropout or z ~ 9 J]^ ]^Q-dropout criteria. 



Searches for star-forming galaxies at high-redshift tend 
to be relatively straightforward in execution. Typically 
these searches involve the acquisition of imaging data 
in three adjacent bands followed by a traditional two- 
color Lyman-Break Galaxy (LBG) selection. While use- 
ful, this approach is limited by the available resources. 
Often times, much larger areas can be searched by utiliz- 
ing suboptimal data sets that still allow for the efficient 
identification of high-redshift star- forming galaxies. 

With the availability of deep wide-area NICMOS i^ieo- 
band data over the two GOODS + HDF-South fields, we 
have such a data set. The wide-area NICMOS + optical 
+ IRAC data allow us to identify sources that have a 
high probability of corresponding to z > 7 galaxies. We 
can perform such a preselection due to the unique colors 
of star-forming galaxies at z > 7. z ^ 7 galaxies are very 
red in the optical to i^ieo colors and moderately blue in 
the HiQo to IRAC colors. 



We then follow up these preselected candidates with 
deep imaging in the Jno-band to better ascertain their 
nature using more traditional two-color LBG criteria. 
For z ^ 7 galaxies (§3.3), we apply the same zs5o — 
Juo I ^110 — ^160 criterion utilized by Bouwens et al. 
(2004), Bouwens & Illingworth (2006), and Bouwens et 
al. (2008). Candidates with very red zgso — ^iio colors 
(corresponding to the Lya cut-off) and blue Juo — i^ieo 
colors (corresponding to the UV continuum) are taken to 
be z ~ 7 galaxies. 

For the z ^7 preselection, we employ Lyman-break- 
like criteria. The candidates must be detected at > 5<j in 
the iJi60 band, have (2^850 — J^16o)ab colors redder than 
1.2, be completely undetected (< 2<j) in the ACS ^435, 
VeoG, and 2775 bands, and not have iJieo — 5.8/im colors 
redder than 2.5 mag (to exclude dust reddened galaxies 
at z~l-2). Our use of a zgso — Him > 1.2 preselection 
was designed to ensure that the candidates we identified 




(^'llO ^16o)aB 

Fig. 2. — 2:850 — ^110 vs. Jiio — i^ieo color-color diagram used 
to select z 7 zsso-dropouts. The orange region shows the lo- 
cation oi z ^ 7 candidates in color-color space. The expected 
colors of high-redshift star-forming galaxies and low-redshift 
galaxy contaminants (shown over the redshift range z ^ 0-2) 
are indicated with blue and red lines, as a function of redshift. 
The plotted uncertainties and limits are la (consistent with 
previous work). The green line enclose the region in color- 
color space where we expect low-mass stars to be found (e.g., 
Knapp et al. 2004) while the red triangles show the position 
of two probable T dwarfs identified within the GOODS fields 
(Bouwens et al. 2008; Mannucci et al. 2006). The cyan points 
show the colors of individual sources in our search fields. The 
lightly-shaded gray region shows the zsso — i^ieo > 1.2 prese- 
lection we use to identify possible z >7 candidates to follow 
up with deeper Jno-band NICMOS observations (§3.2). It is 
apparent that this preselection identifies all eight z ^ 7 can- 
didates (large black squares) previously found by Bouwens et 
al. (2008) and the Bradley et al. (2008) z - 7.6 galaxy. The 
large blue squares show the colors of new sources that satis- 
fied our z ^ 7 zsso-dropout criteria and are included in our 
z ^ 7 Z850-dropout sample (Tables [2] and |4} - while the small 
blue squares show the colors of preselected galaxies that did 
not (and hence are more likely at low redshift). See Table [2] 
for our photometry on those z > 7 candidates identified for 
Jiio-band follow-up observations. 

showed a prominent break across the zgso and i^ieo bands 
(as expected for z > 7 dropouts). However, this criterion 
was not so strong as to exclude any z ^ 7 zgso-dropout 
galaxies identified in previous NICMOS searches (e.g., 
Bouwens et al. 2008; Oesch et al. 2009) or to have a 
substantial impact on the effective selection volume of 
our z ^ 7 search. 

Sources that were detected at l.ba in more than two 
optical bands (^435, Vqog^ or ^775) were also excluded. 
Figure [2] illustrates this preselection in terms of the 
^850 — ^110/^110 — ^160 standard two-color diagram used 
to identify z ^ 7 zgso-dropout galaxies. 

We identified 20 z > 7 candidates for follow-up study 
from the wide- area (~65 arcmin^) NICMOS i^ieo-band 
observations. The coordinates and photometry for candi- 
dates are listed in Tabled After some experimentation, 
we discovered we could fit those z ^7 candidates within 
14 5r X 5r NIC3 fields (i.e., SI, S2, S3, S4, S5, S6, S7, 
S8, Nl, N3, N4, N5, N6, N7). 

We obtained follow-up Jno-band NICMOS observa- 
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Fig. 3. Vboenrs^^sso ^iio^^ieo images (3.5^^ x 3.5^0 of the 7 
z 7 Z850-dropout candidates identified in our new wide-area 
NICMOS data. See Table [2] for the coordinates, magnitudes, col- 
ors, and other properties of the present 2;850 -dropout candidates. 
Deeper Ji25-band observations with the WFC3/IR camera were 
obtained over GNS-zD3 to confirm the weak Jno-band detection 
obtained with NICMOS. The nature of GNS-zD7 is unclear. While 
the source satisfies our zsso-dropout criterion (and would therefore 
appear to be a plausible 2; ~ 7 galaxy), its y — J, J — H colors 
measured with WFC3/IR suggest a red z ^ 1 — 2 galaxy. 



tions on these candidates (typically 3 orbits for each 
NIC3 pointing) with the 60-orbit HST GO program 
11144. Those observations were obtained from Decem- 
ber 19, 2007 to September 9, 2008, with one final or- 
bit of observations taken on October 26, 2009 with the 
WFC3/IR camera. The position of these follow-up data 
is illustrated in Figure [1] with the magenta squares (and 
rectangle for the 1-orbit WFC3/IR observation). 

Once the NICMOS Jno-band observations were ob- 
tained, we redid our photometry on each z >7 candi- 
date. This photometry is included in Table [2] for all 20 
z >7 candidates identified for follow-up study. We then 
applied the 2:350 and Jno-dropout selection criteria we 
describe in the next two subsections to identify probable 
star- forming galaxies at z ~ 7 and z ~ 9, respectively, 
from this preselection. 
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TABLE 3 

Available WFC3/IR photometry for z ~ 7 zsso-dropout galaxy candidates in our wide-area NICMOS sample 



Object ID 


R.A. 


Dec 


^160 


NICMOS 

2^850 — -^110 -^110—^160 


2^850 - ^098 


WFC3/IR 

^098 - -^125 


-^125 — ^160 


GNS-zDl 
GNS-ZD6 
GNS-zD7^ 


03:32:43.29 
03:32:22.66 
03:32:42.84 


-27:42:47.9 
-27:43:00.6 
-27:42:47.7 


25.8±0.2 
25.4±0.2 
26.1±0.3 


1.3±0.4 0.4±0.2 
1.6±0.5 0.2±0.2 
1.5±1.0 0.4±0.3 


0.8±0.4 
1.4±0.4 
0.5±0.9 


0.6±0.3 
0.0±0.2 
0.5±0.6 


0.2±0.1 

0.4±0.1 
0.7±0.3 



^ The measured 2350 — ^105 color ~ 0.0 ±1.0 and J125 — Hieo ~ 0.6 ± 0.2 color suggests this candidate may be a red 2~l-2 galaxy 
and not a 2 ~ 7 star-forming galaxy. 



3.3. z ^ 7 Zs5o-dropout Selection 

We use the same selection criterion for identifying z ~ 
7 Z85o-dropouts as we did in our previous work on the 
HUDF and two GOODS fields (Bouwens et al. 2008). 
Specifically, we require our zgso-dropout candidates to 
satisfy the criteria ((^850 - Jiio)ab) > 0.8) A ((zgso - 
-^iio)ab > 0.8 + 0.4( Jiio - i^ieo)^^) where A represents 
the logical AND symbol. In cases of a non-detection in 
the dropout band, the fiux in the dropout band is set to 
its la upper limit. This two-color selection is illustrated 
in Figure [2] with the position of the z ^ 7 candidates 
from Table [2] included as the large blue squares. 

Z85o-dropout candidates are required to be detected at 
5<j in the i^ieo band (0.6''-diameter aperture) to ensure 
that they correspond to real sources. In addition, in- 
cluded in the present search are also those sources from 
the NICMOS observations from the H. Yan et al. (2010, 
in prep) GO 11192 and Henry et al. (2009) GO 10872 
programs. 

In total, 7 sources from our new search fields satisfied 
our Z85o-dropout criteria. All of these candidates were 
found over the GOODS NICMOS Survey (Conselice et 
al. 2010), and none from the Yan et al. (2010, in prep) 
or Henry et al. (2009) fields. Postage stamps of these 
candidates are shown in Figure [3l Flux measurements 
for the candidates are given in Table [6] in Appendix A. 
The sources range in magnitude from Hiqq^ab ^ 25.4 to 
26.1 mag, with most of the sources being found at ~26 
mag. The surface density of zgso-dropout candidates in 
our fields (~65 arcmin^) brightward of 26.5 mag is ~0.1 
source arcmin"^, very similar to that found by Bouwens 
et al. (2008). 

Overall, the properties ( Jno — ^leo color and sizes) of 
most of our candidates seem consistent with their being 
high-redshift galaxies. The mean half-light radii (includ- 
ing the effect of the PSF) is ~0.25'^ which is similar 
to the candidates found over the HUDF (e.g., Bouwens 
et al. 2008). Meanwhile, the Jno — i^ieo color for the 
sources is 0.6±0.1 mag. This is consistent, albeit some- 
what higher, than the Jno — i^ieo ~ 0.48 color expected 
from our simulations (Bouwens et al. 2008), assuming 
a /3 ~ — 2 (Bouwens et al. 2010a; Stanway et al. 2005; 
Finkelstein et al. 2010; Bunker et al. 2010). Nonetheless, 
a few of the candidates (GNS-zD3, GNS-zD5) in our se- 
lection have sufficiently red Jno — i^ieo colors (~0.8 to 
1.0 mag) that they might be low-redshift interlopers. 

The NICMOS Jno-band observations available on 
GNS-zD3 were only moderately deep (~2 orbits) and 
so GNS-zD3 only shows a weak detection in that band. 
Fortunately, we were able to obtain deeper Ji25-band ob- 
servations on GNS-zD3 with the WFC3/IR instrument 



V i z J H 




Fig. 4. V^oenrs 2^850 -^iio^ieo images (3.5^^ x 3.5^0 of the 2 
sources in our new wide-area NICMOS data that satisfy our z 9 
Jiio-dropout criteria. However, we consider it unlikely that ei- 
ther of the candidates identified here corresponds to a zr^9 galaxy. 
GNS-JDl appears to be detected at la in both the Veoe and 2775 
bands, not sufficient for us to rule it out as a 2; ~ 9 Jno-dropout 
candidate but suggesting it may be a 2;~l-2 galaxy. GNS-JD2, 
by contrast, shows no evidence for being detected at wavelengths 
other than 1.6/im (if-band). This may suggest it is a transient 
source (SNe: see §3.4) or spurious (since it is close to the edge of 
the NIC3 field where it was found). However, we cannot rule out 
the possibility it corresponds to a 2; ~ 9 galaxy (but we consider it 
very unlikely). See Tables [2] and [6] for the coordinates, magnitudes, 
colors, and other properties of these candidates. 

(§2). The new data confirms there is a detection in the 
J band, and the measured magnitude 26. 5 ±0.2 is consis- 
tent with what is found with NICMOS. 

Three of our new zgso-dropout candidates (GNS-zDl, 
GNS-zD6, GNS-zD7) are found over the upper region in 
the CDF-South where deep wide-area WFC3/IR obser- 
vations were recently taken as part of the Early Release 
Science program (G011359: PI O'Connell). Given that 
these observations reach ~0.7-1.0 mag deeper than the 
NICMOS observations utilized in this study and extend 
over three bands loos, ^125, and i^ieo, they are useful 
for characterizing the typical dropout candidates found 
in this search. 

What do these deeper data suggest about the candi- 
dates in our selection? Photometry on the three afore- 
mentioned Z85o-dropout candidates was performed using 
the new WFC3/IR observations (utilizing specifically the 
Bouwens et al. (2010c) reductions). The results are sum- 
marized in Table m and the conclusions are mixed. GNS- 
zD6 is clearly a 2: ~ 6 galaxy - though from the measured 
colors its redshift is likely in the range z~6.2-6.5. The 
nature of GNS-zDl is slightly less clear from the data. 
Its measured z — J — J colors are consistent with its 
being either a z ~ 6.5 galaxy or a red z^l-2 galaxy. The 
blue J - H - 0.1 colors of GNS-zDl seem to slightly 
favor the case that it is a z > 5 galaxy. Finally, GNS- 
zD7 seems most consistent with being an intrinsically red 
z^l-2 galaxy, having very red J — H ^ 0.6, Y — J ^ 0.7 
colors. 
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These results suggest that the present zsso-dropout se- 
lection is successful in identifying z > 6.5 galaxies, albeit 
with a mean redshift somewhat lower than the z ~ 7.3 
estimated in Bouwens et al. (2008). The lower mean 
redshift for the sample is consistent with the expected 
bias based upon evolution across the zgso-dropout selec- 
tion window (where more luminous galaxies are present 
e.g. at z ~ 6.5 than at z ~ 8: Muhoz & Loeb 2008). 
The likely contamination of our selection by one proba- 
ble low-redshift source (GNS-zD7) is consistent with the 
24% contamination rates estimated in §3.6. 

The principal reason we are finding modest levels of the 
contamination over the GOODS fields is because of the 
limited depth of the available ACS optical data over the 
GOODS fields. This contamination is somewhat higher 
than estimated over other z ^ 1 zgso-dropout selections 
like the HUDF - where it was estimated to be ~12% (e.g., 
Bouwens et al. 2008). In order to reduce these contami- 
nation levels, it would therefore be ideal if deeper optical 
data - particularly in the F606W, F775W, and F814W 
bands - could be obtained over the GOODS fields. 

3.4. z ~ 9 Jiio-Dropout Selection 

Similar to our z ^ 7 zgso-dropout selection, we adopt 
the same z ^ 9 Jno-dropout selection criteria as we 
used in the -23 arcmin^ Bouwens et al. (2008) NICMOS 
search. Jno-dropout candidates in our selection are re- 
quired to satisfy the criterion ( Jno — Hiqo)ab > 1.3 and 
not show > 2(7 detections in any of the optical bands 
(or > 1.5(7 in two bands), z ^ 9 Jno-dropout candi- 
dates are also required to be detected at 6a in the i^ieo 
band (0.6''-diameter aperture) to ensure that most of the 
sources are real. We use a 6(7 detection criterion for our 
Jiio-dropout selection (instead of a 5(7 criterion) because 
we only have one passband to evaluate the reality of the 
candidates. 

In total, we identified 2 sources that satisfied our z ^ 9 
Jiio-dropout criteria. Postage stamps of these candi- 
dates are provided in Figure HI and their photometry is 
summarized in Table [2l Flux measurements for the can- 
didates are given in Table [6] in Appendix A. 

While both candidates formally satisfy our Jno- 
dropout selection criteria and therefore may correspond 
to z^9 galaxies, neither candidate is very compelling, 
and there are reasons to suspect each candidate may be 
a contaminant. For example, GNS-JDl is formally de- 
tected at 1(7 in both the Veoe and ^775 bands - suggesting 
that it may actually correspond to a 2: — 2 Balmer-break 
galaxy. In addition, GNS-JDl is detected at 5(7 in both 
the 3.6/im and 4.5/im bands, with ma.e/im = 25.1 ± 0.2 
mag and m^^^j^m = 25.0 ± 0.2 mag. This latter IRAC 
photometry was performed by modelling the light pro- 
files of nearby neighbors in the IRAC imaging data, sub- 
tracting these model profiles from the observations, and 
then measuring the fiux in simple 2.5''-diameter aper- 
tures (e.g., Labbe et al. 2006; Gonzalez et al. 2010). 

The other Jno-dropout candidate GNS-JD2 is possibly 
spurious. None of the observations at other wavelengths 
provide any evidence it is real. It is not detected in the 
ACS observations, the NICMOS Jno-band observations, 
nor even the IRAC 3.6/i and 4.5/i observations. One 
would have expected GNS-JD2 to be formally detected 
in the IRAC 3.6/i data at ~1. 5-2(7 if it was well sepa- 
rated from its neighbors. However, since the source is not 



well separated, fiux measurements are sufficiently chal- 
lenging that the IRAC non-detection does not definitely 
argue against the reality of the source. Nonetheless, the 
non-detection does suggest that the original detection of 
GNS-JD2 in the i^ieo-band may have been anomalous 
and therefore the source is likely either spurious or a 
transient source (i.e., a SN). 

Additional evidence for the source being spurious 
comes from its proximity to the edge of the NIC3 field 
where it was found (due to enhanced non-Gaussianity of 
the noise there). To test the spurious hypothesis, we ex- 
amined the individual NICMOS exposures that went into 
its i^i6o-band stack. None of these exposures provided a 
dominant contribution to the cumulative i^ieo-band fiux, 
leaving us with no compelling reason to fiag this particu- 
lar candidate as spurious. An alternate hypothesis is that 
GNS-JD2 is a SN. As we discuss in §3.6, we would expect 
~1 SN to be identified over our ~65 arcmin^ i^ieo-band 
observations as a z > 7 galaxy candidate. 

Besides these 2 formal Jno-dropout candidates, there 
were 3 other z > 7 candidates identified in the NIC- 
MOS i^i6o-band data for follow-up (GNS-Spl, GNS- 
Sp2, GNS-Sp3 in Table [2]), but which were not detected 
in the later Jno-band observations. This would give 
these sources nominal Jno — i^ieo colors > 2 mag (Icr) 
and make them possible z ^ 9 Jno-dropout candidates. 
However, there are reasons to be concerned about these 
candidates and whether they correspond to real sources. 
While the formal significance of each source is ~5(7 in 
the i^i6o-band data, these sources are detected in stacks 
of <9 NICMOS exposures (and hence subject to noise 
with significant non- Gaussian signatures) and are typi- 
cally near a bright source (GNS-Spl, GNS-Sp3). Con- 
sequently, the true significance of these sources is some- 
what smaller, i.e., ^4(7, making it more likely they are 
spurious (see also the simulation results in §3.6) O These 
three sources are included in Table [2] under the label 
"Probable Spurious J-dropout Candidates." 

3.5. Total Sample of z ^ 7 Zs5o-dropouts and z ^ 9 
Jiio- dropouts 

In §3.3 and §3.4, we identify 7 plausible z ^ 7 z^bq- 
dropout candidates in the deep, wide- area (~65 arcmin^) 
NICMOS i^i6o-band data recently obtained over the 
GOODS + HDF-South fields. Two possible z - 9 Jno- 
dropout candidates are identified, but appear unlikely 
to correspond to z^9 galaxies. To increase our search 
area to ~ 88 arcmin^ (including all the NICMOS ob- 
servations in Table [1]), we combine this sample with 
the sample of z > 7 dropouts found in the deep, but 
smaller area ~23 arcmin^ NICMOS data already consid- 
ered by Bouwens et al. (2008). Figure [5] shows postage 
stamp images of the z ^ 7 field candidates from that 
study. Those postage stamps are essentially identical to 
those presented in Bouwens et al. (2008), but incorpo- 
rate 8 additional orbits of NICMOS observations taken 
on two HUDF ^850-dropout candidates (UDF-387-1127, 

Note that this is in contrast to situations where similar signifi- 
cance sources are found in observations created from a much larger 
number of exposures (e.g., the 56 exposures stacks used for the 
WFC3/IR HUDF09 ifieo-band observations: Oesch et al. 2010a; 
Bouwens et al. 2010b). In those cases, the noise characteristics 
are much closer to Gaussian, and apparent 5cr sources are indeed 
significant at the Scr level. 
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TABLE 4 

Our total sample of 2; ~ 7 2;85o-dropout galaxy candidates selected from NICMOS observations over the CDF-S + HDF-N 

GOODS FIELDS, THE HUDF, AND GALAXY CLUSTER FIELDS (SEE TABLE (TJ.^ 



Object ID 


R.A. 


Dec 


^160 


2^850 — JllO 


Jiio — HiQo 


ifl60 — Ks 


Muv,AB^ 


Ref^ 




Candidates in the CDF-South 


and HDF-North Fields 














12:36:44.68 


62:16:15.4 


25.4±0.2 


1.7±0.8 


1.0±0.3 




—21.5 








GNS-zD6 


03:32:22.66 


-27:43:00.6 


25.4±0.2 


1.6±0.5 


0.2±0.2 




-21.5 


[12,131 


GNS-ZD4 


12:36:10.93 


62:09:15.6 


25.8±0.3 


>1.6 


0.7±0.3 




-21.1 








GNS-zDl 


03:32:43.29 


-27:42:47.9 


25.8±0.2 


1.3±0.4 


0.4±0.2 





-21.1 








HDFN-3654-1216 


12:36:54.12 


62:12:16.2 


26.0±0.1 


1.1±0.3 


0.4±0.3 


0.0±0.3 


-20.9 


[2,71 


GNS-ZD3 


03:32:06.10 


-27:46:37.3 


26.1±0.3 


1.5±1.0 


0.8±0.5 


— 


-20.8 








GNS-zDT^ 


03:32:42.84 


-27:42:47.7 


26.1±0.2 


1.5±1.0 


0.4±0.5 





-20.8 








UDF-640-1417 


03:32:42.56 


-27:46:56.6 


26.2±0.1 


1.3±0.3 


0.5±0.2 


0.5±0.3 


-20.7 


[1,2,3,4, 


5, 


7,11,12] 


GNS-ZD2 


03:32:32.03 


-27:45:37.2 


26.2±0.3 


1.6±1.0 


0.6±0.3 




-20.7 








NICPAR2-3308-5229^ 


03:33:08.29 


-27:52:29.2 


26.7±0.1 


> 2.2^ 


0.6±0.2 




-20.2 


[7,8] 


CDFS-3225-4627 


03:32:25.22 


-27:46:26.7 


26.7±0.2 


1.4±0.5 


0.2±0.3 


<0.2f 


-20.2 


[2,7] 


UDF-983-964 


03:32:38.80 


-27:47:07.2 


26.9±0.2 


>3.2^ 


0.1±0.3 


0.0±0.8 


-20.0 


[1,2,4, 


5, 


7,8,11] 


UDF-387-1125S 


03:32:42.56 


-27:47:31.4 


27.1±0.2 


1.3±0.5 


0.8±0.4 


<o.of 


-19.8 


[1,5 


,7 


,11] 


UDF-3244-4727g 


03:32:44.02 


-27:47:27.3 


27.3±0.2 


>2.6f 


0.5±0.4 


<o.if 


-19.6 


[7,^ 




NICPARl-3303-4111^ 


03:33:03.81 


-27:41:12.1 


27.8±0.1 


>1.5f 


0.4±0.2 




-19.1 












Candidates 


in Lensing 


Cluster Fields'" 














A1689-ZD1 


13:11:29.73 


-01:19:20.9 


24.7±0.1 


>2.2f 


0.6±0.2 




-19.8 




6 




CL0024-ZD1 


00:26:37.93 


17:10:39.0 


25.6±0.1 


1.3±0.8 


0.4±0.3 




-19.3 
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CL0024-iDl 


00:26:37.78 


17:10:40.0 


25.0±0.1 


0.9±0.2 


O.liO.l 




-19.9 
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^ There are a few weaker 2850-dropout candidates behind lensing clusters that have been identified by Richard et al. (2008) and Bouwens et al. 
(2009a), but none of these have deep enough optical data to be included in the present list as reliable 2 ~ 7 candidates. 

^ The absolute magnitudes estimated for the sources are for an effective rest-frame wavelength of ~1900A. The absolute magnitudes given for the 
lensed sources include a correction for the estimated magnification factor (Bradley et al. 2008; Zheng et al. 2009). 

References: [1] Bouwens et al. (2004b), [2] Bouwens & Illingworth (2006), [3] Yan & Windhorst (2004), [4] Coe et al. (2006), [5] Labbe et al. 
(2006), [6] Bradley et al. (2008), [7] Bouwens et al. (2008), [8] Oesch et al. (2009), [9] Zheng et al. (2009), [10] Gonzalez et al. (2010), [11] Oesch 
et al. (2010a), McLure et al. (2010), Bunker et al. (2010), Yan et al. (2010), Finkelstein et al. (2010), [12] Castellano et al. (2010), [13] Hickey et 
al. (2010) 

^ This source has a measured zgbo — ^105 color ~ 0.0 ± 1.0 and J125 — Hiqo ~ 0.6 ± 0.2 color in the high S/N WFC3/IR observations (see Table [3|) . 
This suggests this candidate may be a red 2~l-2 galaxy and not a 2 ~ 7 star-forming galaxy. 
^ See the footnotes in Table 2 of Bouwens et al. (2008) for these candidates 
^ Upper and lower limits on the measured colors are the la limits. 

S The S/N on our Jno-band fiuxes for UDF-387-1127 and UDF-3244-4727 is higher than initially reported in Bouwens et al. (2008). The improved 
S/N is the result of 8 orbits of additional Jno-band observations on these sources. 

UDF3244-4727) in the Jno band. 

Our total sample of z ^ 7 zgso-dropouts identified in 
our search fields (Table [T]) - including both the old and 
new data - is presented in Table IH Also included in this 
table are three z ^ 7 zgso-dropouts identified in searches 
behind lensing clusters (Bradley et al. 2008; Bouwens et 
al. 2009; Zheng et al. 2009). 



3.6. Estimated Contamination Rate 

High redshift dropout selections are subject to contam- 
ination from (1) sources that enter the selection due to 
photometric scatter, (2) transient sources, (3) low-mass 
stars, and (4) spurious sources (i.e., corresponding to no 
real source). We consider each of these sources of con- 
tamination in examining the present z ^ 7 zgso-dropout 
and z ^ 9 Jno-dropout selections. 

Possible Contamination from Photometric Scatter: In 
general, the most important contaminant for high- 
redshift dropout selections are sources that enter the se- 
lections due to photometric scatter (e.g., see Bouwens 
et al. 2008; Bouwens et al. 2010b). We estimate the 
contamination rate that results from this effect by start- 
ing with the color distribution observed for i^~25.0-26.5 
sources in the HUDF09 WFC3/IR field (Bouwens et al. 
2010b; Oesch et al. 2010a) and then adding photomet- 
ric scatter to match that expected for each source in our 
NICMOS catalogs. Repeating these simulations for each 
source in our NICMOS catalogs (from all search fields 



not considered in Bouwens et al. 2008), we predict that 
1.7 sources and 0.2 sources would enter our zgso-dropout 
and Jiio-dropout selections, respectively, via photomet- 
ric scatter. The implied contamination rate is 24% and 
10%, respectively, which is somewhat higher than esti- 
mates for our HUDF 2:85o-dropout selections (Bouwens 
et al. 2008) or for our new WFC3/IR results (Oesch et 
al. 2010a; Bouwens et al. 2010b). This is due to the 
somewhat shallower depths of the ACS optical data over 
the GOODS fields and larger prevalence of low-redshift 
sources with similar colors to high-redshift galaxies at 
bright magnitudes (see also §3.3). The above procedure 
is essentially identical to that used in Bouwens et al. 
(2008). 

Possible Contamination from Transient Sources: Con- 
tamination from transient sources, particularly SNe, are 
potentially important for each of our search fields, given 
that the NICMOS observations were typically acquired 
at least 4 years after the ACS optical observations. Us- 
ing the GOODS SNe searches (Riess et al. 2004; Strolger 
et al. 2004) as a baseline, Bouwens et al. (2008) argued 
that the contamination rate from SNe should be no larger 
than 0.012 arcmin"^ for observations where the optical 
and near-IR J -\- H observations are taken at very dif- 
ferent times. This suggests a contamination rate of ~1 
source for the present search for z > 7 galaxies over ~65 
arcmin^ of new data. However, the above calculation as- 
sumes the J and H band observations are acquired at the 
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Fig. 5. — Veoenrs^^ssO'/iio^ieo images of the 8 2; ~ 7 2:850- 
dropout candidates previously identified in the ultra-deep, wide- 
area NICMOS data (Bouwens et al. 2008: but see also e.g. Oesch 
et al. 2009), but now utilizing the deeper Jno-band observations (8 
additional orbits) obtained on two candidates (UDF-3244-4727 and 
UDF-387-1125) as part of the GO 11144 program. The Jno-band 
detections for these two candidates are now much more significant 
than they were in Bouwens et al. (2008: Figure 3 from that work). 
The other 6 candidates are included here for completeness. Each 
of the above candidates is detected at > 4.5cr in both the Jno and 
ifi60 bands. 

same time. In reality, for most of the new data, the Jno- 
band observations were taken at least 3-12 months af- 
ter the NICMOS i^ieo-band observations for each of our 
new Z85o-dropout candidates. Consequently, the mea- 
sured J-H colors for SNe's found in GOODS NICMOS 
survey would likely be very red and hence show up as an 
apparent Jno-dropout candidate. We might therefore 
expect 1 SNe contaminants to be present in our Jno- 
dropout selection. This may be the case for GNS-JD2 
(see §3.4). 

Possible Contamination from Low-mass Stars: Low- 
mass stars (e.g., T or L dwarfs) have very similar z — J, 
J — H colors to star- forming galaxies at z ^ 7 and there- 
fore could contaminate our zgso-dropout selection. Given 
the modest resolution of the NICMOS data, it is difficult 



to use these data to clearly distinguish small (~0.15'' 
half-light radius) star-forming galaxies found at z ~ 7 
from unresolved low-mass stars. Fortunately, the ^850- 
band observations have sufficient resolution and depth 
that low-mass stars (with z — J colors ~1.5-2) would be 
evident at 5a (0.2''-diameter aperture) in the zgso-band 
images as point sources. Only one possible T dwarf can- 
didates is evident in the new wide-area NICMOS obser- 
vations and that is the source at 03:32:22.66, —27:43:00.6 
(GNS-zD6). However, that source appears to be resolved 
in the new WFC3/IR observations as part of the Early 
Release Science Program (G011359: PI O'Connell), so 
that source appears unlikely to be a low mass star. The 
present situation is somewhat in contrast to the NIC- 
MOS observations considered by Bouwens et al. (2008) 
where two probable T-dwarfs were identified in selecting 
z ^ 7 galaxies (included on Figure [2] as the red triangles). 

Possible Contamination from Spurious Sources: Con- 
tamination from spurious sources could be a concern for 
our z ^ 9 Jno-dropout selection. Each of our Jno- 
dropout candidates is only detected in a single band and 
the estimated ~6(T significance of these detections may 
be an overestimate (due to real data possessing many 
non-Gaussian characteristics). Therefore, to estimate 
the likely number of spurious sources, we used the stan- 
dard negative image test (e.g., Dickinson et al. 2004) and 
repeated our z ^ 9 Jno-dropout selection on the negative 
iJi6o-band images (after masking out sources brighter 
than 23 mag). Such a selection yielded ~12 Ji25-dropout 
candidates, but all of the candidates (on the negative im- 
ages) could be eliminated due to an obvious association 
with defects or irregularities in the reduction. While this 
might suggest that spurious sources do not dominate our 
Jno-dropout selection, there were a modest number of 
formal 5cr detections on the negative images quite close 
to satisfying our Jno-dropout selection criteria, and so 
contamination from spurious sources here is certainly not 
impossible. 

By contrast, for our z ^ 7 zgso-dropout selection, spu- 
rious sources are not an important concern. Each of our 
candidates is detected at >2a in the Jno band and >5a 
in the iJieo band - making contamination from spurious 
sources extremely unlikely. 

Summary: In total, we expect 1.7 and 1.2 contami- 
nants in our z ^ 7 zgso-dropout and z ^ 9 Jno-dropout 
selections, respectively (equivalent to contamination lev- 
els of 24% and 60%). The only meaningful source of 
contamination for our zgso-dropout selections is photo- 
metric scatter, while for our Jno-dropout selection, sev- 
eral sources of contamination contribute. We expect ~1 
contaminant from transient sources (SNe), 0.2 contam- 
inants from photometric scatter, and <1 contaminants 
from spurious sources. 



4. 



IMPLIED CONSTRAINTS ON THE REST-FRAME UV 
LUMINOSITY FUNCTION 



The present wide- area (~88 arcmin^) NICMOS search 
for z > 7 candidates provides a useful constraint on 
the volume density of luminous star-forming galaxies at 
z ^ 7 and z ^ 9. These observations cover twice as much 
area as ~39 arcmin^ ERS observations and cover simi- 
lar area to that available from the CDF- South HAWK-I 
observations (90 arcmin^: Castellano et al. 2010; Hickey 
et al. 2010). They therefore provide a somewhat compa- 
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(left) Present determination of the stepwise UV LF at z ~ 7 using deep wide-area NICMOS + ground-based observations (large 



solid red circles: see §4). The z 7 UV LF is only redetermined here at luminosities where we can take advantage of the new ~65 arcmin^ 
of NICMOS observations (see Table [T} to search for candidate z >7 galaxies. For comparison, we also show the z 7 LFs reported by 
Bouwens et al. (2008: hlack crosses), Ouchi et al. (2009: black squares), Castellano et al. (2009: green squares), Oesch et al. (2010a: open 
magenta circles), McLure et al. (2010: blue squares), and Bouwens et al. (2010c: open red circles). Constraints from Oesch et al. (2009) 
and Wilkins et al. (2010b) are similar, but not shown to reduce confusion. The solid red line is the best-fit Schechter function presented 
in Oesch et al. (2010a). (right) Stepwise UV LF determined here at z ~ 7 (solid red circles) versus that derived at z ~ 4 (blue), z ~ 5 
(green), and z ^ 6 (cyan) by Bouwens et al. (2007). The open red circles and red line represent the z ^ 7 LF determined by Bouwens et al. 
(2010c). Also shown are our constraints on the LF at z^9 from the present Jno-dropout search (black upper limits: see §4). Similar to the 
situation for our z ~ 7 LF, the z ^ 9 Jno-dropout LF is only computed at luminosities where we can take advantage of the new NICMOS 

^L*^3 (~— 21 mag) star-forming galaxies at z ~ 7 and z ~ 9 is IS^^x and >25x lower, respectively, than 



data. The volume density of 
at z ~ 4. 



rable constraint on the bright end (i.e., Muv < —20.7) 
of the z > 7 LFs to these studies. And while the Ouchi 
et al. (2009) Subaru Suprime-Cam search extends over 
considerably more area (1568 arcmin^) than these obser- 
vations, the present NICMOS search is somewhat deeper 
than that study (by ^0.7 mag). The present search also 
benefits from somewhat higher quality multiwavelength 
data - deep high-resolution ACS or mid-IR IRAC data 
(Giavahsco et al. 2004; Dickinson & GOODS Team 2004) 
- than are generally available for the Subaru fields. This 
allows for a more robust discrimination between high- 
redshift star-forming galaxies and low-redshift contami- 
nants (and low mass stars). 

Here we will estimate the stepwise LF <I>(M) using 
the relatively simple V/Veff approach (e.g., Steidel et 
al. 1999) where 



<l>(m) = 



N{m) 

Veff{m) 



(1) 



where N{m) is the number of z > 7 candidates in a given 
magnitude interval m (from all of our search fields) and 
Veff{m) is the effective volume in the magnitude interval 
m. 

To compute the effective volume with which we can se- 
lect star- forming galaxies in our various search fields, as 
a function of magnitude and redshift, we run detailed 
Monte- Carlo simulations where we introduce model 
galaxies into the observed data and attempt to select 
them as dropouts using the procedures laid out in §3. 
These simulations are performed in the same way as they 
were performed in Bouwens et al. (2008). For the model 
galaxies in these simulations, we assume that the z > 7 
galaxies of a given luminosity have similar pixel-by-pixel 
morphologies to z ^ 4 galaxies in the HUDF with the 



TABLE 5 
Stepwise Constraints on the 
rest-frame uv lf at z ~ 7 and z ~ 9 
FROM Wide- Area NICMOS + 
Ground-Based Observations (§4).^ 

Muv,AB^ cl>k (10-^ Mpc-3 mag-i) 

zsso-dropouts (z ~ 7) 
-21.55 8^5^ 
-20.95 86^^^ 
J-dropouts (z ~ 9) 
-21.85 < 11^ 

-21.25 < 22^ 

^ The UV LFs are only redetermined here 
at luminosities where we can take advan- 
tage of the new ~65 arcmin^ of NICMOS 
observations (see Table [T} to search for 
candidate z >7 galaxies. 
^ The effective rest-frame wavelength is 
~1900A for our zsso-dropout selection 
and ~1500A for our Jno-dropout selec- 
tion. 

^ Upper limits here are la (68% confi- 
dence). 

same luminosity, but scaled in size as (1 + 2:)"^ to match 
the observed size-redshift trends at 2: > 2 (e.g., Fergu- 
son et al. 2004; Bouwens et al. 2004; Bouwens et al. 
2006; Oesch et al. 2010b). The /7V-continuum slopes 
of the model galaxies are assumed to have a mean of —2, 
with a la scatter of 0.5. These latter assumptions match 
the colors of bright star-forming galaxies found at 2: ~ 7 
(Bouwens et al. 2008; Oesch et al. 2010a; Bunker et al. 
2010; Bouwens et al. 2010b; Finkelstein et al. 2010). 

In computing UV LFs on the basis of our 2: > 7 search, 
we consider all seven z ^ 7 2;85o-dropout candidates in 
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the newest NICMOS data and suppose that ~5.3 2^850- 
dropout candidates from this search correspond to 2: ~ 7 
galaxies (accounting for the 24% contamination rate es- 
timated in §3.6). We also incorporate the constraints 
from the NICMOS + ground-based z ^ 1 search consid- 
ered by Bouwens et al. (2008). In addition, we suppose 
that <1 Jiio-dropout candidates from this search cor- 
respond to galaxies, given the concerns that exist 
for each candidate (§3.4). Finally, the z ^ 1 and 2: ~ 9 
LFs are only redetermined here at luminosities where we 
can take advantage of the new ~65 arcmin^ of NICMOS 
observations (see Table [I]). 

Our stepwise UV LF at 2: ~ 7 is presented in Table [5] 
and in Figure [6] (/e/i^). For comparison. Figure [6] also in- 
cludes the z ^ 1 LFs of Bouwens et al. (2008), Ouchi et 
al. (2009), Castellano et al. (2010), Oesch et al. (2010a), 
McLure et al. (2010), and Bouwens et al. (2010c). The 
present LF results are in reasonable agreement with pre- 
vious determinations. In the right panel of Figure [6l the 
present z ^ 1 LF results are shown relative to the LFs at 
z^A^z^b^z^Q (from Bouwens et al. 2007) to provide 
a sense of the evolution from 2; ~ 4. Also included on 
this figure {black upper limits) are the constraints on the 
UV LF at 2; ~ 9 from the present Jno-dropout search. 
The UV LF is 13^5 x lower at ^ - 7 than at ^ - 4 and 
>25x lower (la) at 2: ~ 9 than at 2: ~ 4. The latter 
constraint is the most stringent constraint yet available 
on the volume density of > L*=2, 9(^^(^x'i^s at z ^ 9. 

Of course, the above LF determinations are subject to 
uncertainties as a result of large-scale structure ( "cosmic 
variance"), and therefore to properly frame the above 
results, it is helpful to estimate the size of these uncer- 
tainties. For convenience, we utilize the Trenti & Sti- 
avelli (2008) cosmic variance calculator to make this esti- 
mate. Given that the approximate volume density of the 
sources probed at z ~ 7 and z ~ 9 is -6x10"^ Mpc-^ 
and ~2xl0~^ Mpc~^, respectively, the approximate bias 
parameters are 10 and 12, respectively (e.g., Somerville 
et al. 2004; Trenti & Stiavelh et al. 2008). Given that 
our NICMOS search data are fairly randomly scattered 
over the two GOODS fields, we assume two independent 
search fields of dimension 10' x 16' and adopt a width for 
our redshift selection window Az ~ 1.5 (e.g., see Figure 
7 of Bouwens et al. 2008). Utilizing these inputs, we es- 
timate that the large-scale structure uncertainties on our 
LF results here are ~17% and ~20%, respectively. Large- 
scale structure uncertainties are therefore much smaller 
than those uncertainties estimated from the small num- 
ber statistics (^8 sources) and uncertain contamination 
rates (§3.6). 
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5. SUMMARY 

We have taken advantage of ~88 arcmin^ of deep, wide- 
area NICMOS data to search for 2: > 7 galaxies within 
the HUDF, the two GOODS fields, and the HDF South. 
This search incorporates ~65 arcmin^ of wide-area NIC- 
MOS data not previously used to identify z > 7 galax- 
ies0 ~248 arcmin^ of deep ground-based data (> 25.5 
mag, 5<j) previously used by Bouwens et al. (2008) is also 
considered. In total, we find ~7 plausible z ^ 7 zg^o- 
dropout candidates in the new NICMOS observations 
(six of which are being reported for the first time) and 
~2 possible (but probably unlikely) z ^ 9 Jno-dropout 
candidates. These candidates significantly add to the 
number of luminous z ^ 7 candidates known within the 
GOODS fields and improve our constraints on the volume 
density of luminous galaxies at z ^ 7. These candidates 
have recently been used to model the stellar populations 
of bright z ^ 7 galaxies (Gonzalez et al. 2010; Labbe et 
al. 2010b). 

When taken together with the NICMOS data already 
considered in Bouwens et al. (2008), we have identified 15 
z ^ 7 2:85o-dropout candidates in total from ~88 arcmin^ 
of deep NICMOS data. After running detailed simula- 
tions to estimate the selection volumes and contamina- 
tion rates, we use our new z ^ 7 samples (plus ground- 
based search area) to update the Bouwens et al. (2008) 
determination of the UV LF at z ~ 7 and to strengthen 
our constraints on the LF at z ^ 9. 

We find that the bright end of the UV LF at z ~ 
7 is 13^5 X lower at z ^ 7 than at z ~ 4. At z ^ 
9, the UV LF is a factor of >25x lower (la) than at 
z ~ 4, assuming that at most one of the Jno-dropout 
candidates identified here is at 2: ~ 9. 
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TABLE 6 

Measurements of the fluxes of our z>7 candidates in the optical 

ACS AND NEAR-IR NICMOS DATA.^ 



fi^ (nJy) 



Object ID 
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-^110 
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APPENDIX 

FLUX MEASUREMENTS FOR OUR INDIVIDUAL Z>7 CANDIDATES 

To better assess the significance of the colors measured for our z^7 candidates, we also tabulate the observed fluxes 
of each candidate in Table [6l The tabulated fluxes also allow us to readily evaluate whether any of the candidates is 
detected in the optical (and therefore not likely at high redshift). 



